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Abstract 

In i? 2 -inflation scalaron slow roll is responsible for the inflationary stage, while 
its oscillations reheat the Universe. We find that the same scalaron decays induced 
by gravity can also provide the dark matter production and leptogenesis. With R 2 - 
term and three Majorana fermions added to the Standard Model, we arrive at the 
phenomenologically complete theory capable of simultaneously explaining neutrino os- 
cillations, inflation, reheating, dark matter and baryon asymmetry of the Universe. 
Besides the seesaw mechanism in neutrino sector, we use only gravity, which solves all 
the problems by exploiting scalaron. 

1 Introduction and summary 

The first inflationary model widely discussed in literature, dubbed i? 2 -inflation [1], works 
in a very economic way, exploiting one and the same interaction — gravity — to accomplish 
both inflation and subsequent reheating. It is tempting to exploit gravity somewhat further, 
addressing two other important issues: dark matter production and generation of the baryon 
asymmetry of the Universe. 
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To this end we consider in this work the universal mechanism of particle production 
operating in i? 2 -inflationary model: scalaron decay. We have found that free scalars heavier 
than 10 keV and fermions heavier than 10 7 GeV are forbidden in this scenario, as they would 
overdose the Universe. The fermion of 10 7 GeV is a viable dark matter candidate in this 
model, while light scalars could contribute to the hot dark matter component at best. 

We have also found that with two additional right-handed sterile neutrinos one can 
explain the neutrino oscillations (via standard seesaw mechanism) and baryon asymmetry 
of the Universe (via standard non-thermal leptogenesis) . Curiously, the amount of baryon 
asymmetry available in the model is strongly constrained from above, so that the observed 
amount is only one order of magnitude below the model upper limit. 

Both the dark matter fermions and sterile neutrinos are produced in post-inflationary 
Universe in scalaron decays. This suggests that SM, supplemented with i? 2 -term and three 
right-handed sterile neutrinos, forms a kind of naturally complete theory. This theory ex- 
plains different phenomena beyond the SM — inflation, reheating, dark matter, baryon asym- 
metry of the Universe — involving one and the same mechanism based on the peculiarities of 
scalaron interactions with itself and other fields. These phenomena together with neutrino 
oscillations are the main observational facts pointing at incompleteness of the SM. All of 
them can be explained within the proposed model. 

2 Gravitational production of Dark Matter in i? 2 -inflation 

We start with the following Lagrangian in the Jordan frame: 



where we use the reduced Planck mass M P related to the Planck mass M Pi as M P = 
Mpij\f%T\ = 2.4 x 10 18 GeV; S^ atter includes the action of the Standard Model and other 
new fields. In particular, for free scalar <p and Dirac fermion ip one has 1 



1 These fields are free in the sense of particle physics: in both the Jordan (Eqs. (2) and (3)) and Einstein 
frames (Eqs. (5) and (6)), they are free at Mp — > oo. 




(1) 





(3) 
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where T> is the covariantly generalized Dirac operator, see e.g. [4]. It is convenient to go to 
the Einstein frame by the conformal transformation 

9vw -> 9vu = XQ^u , X = exp <f>/Mp) . 

Scalar and fermion fields are rescaled then as 

<p = x ~ 1/2 <f, 4> -+ i= = x" 3/4 ^ , v v = x ~ 1/2 v ■ 

In the Einstein frame the gravity action takes the Einstein-Hilbert form, but additional 
scalar degree of freedom emerges and couples to all matter fields. Thus the original action 
(1) transforms into 



iEF 



= \/^gd A x R +- g^da <pd u (p — - 1 — 



+ ^matter i (^) 



where S^ tter includes interactions with the field 0. In particular, now the actions (2) and 
(3) read 

'(5) 

s$ F = Jy/H d'x (i$b* - ^| ■ (6) 

At small values of this new degree of freedom decouples from all other fields and both 
frames become identical, as in this limit R is also small. 

Cosmology of the homogeneous and isotropic Universe described by the action (1) has 
inflationary stage [1] at large values of R. In the Einstein frame this stage is realized as 
large-field inflation in the slow-roll regime taking place at super-Planckian values of the field 
which serves as the inflaton. The equivalent scalar mode in the action (1) was named 
scalaron [1] and we use both names in what follows. 

Inflation gives rise to primordial scalar perturbations, whose amplitude normalization to 
the observed CMB anisotropy and large-scale structure yields the estimate [2] 

fi = 1.3 x 1(T 5 M P . 

The spectral index of the scalar perturbations and parameters of the generated tensor per- 
turbations are consistent with observational constraints [3]. 
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When the slow-roll conditions get violated, inflation terminates and the inflaton starts 
to oscillate rapidly, with frequency equal to scalaron mass This drives the Universe expan- 
sion like at matter-dominated stage. The intermediate stage naturally ends up with inflaton 
decays into ordinary particles due to universal interactions as in (5), (6). In particular, 
scalaron decay rates into a pair of sufficiently light scalars and into a pair of sufficiently light 
fermions are 2 

= 192^ M 2 ' (7) 
_ _ ym% 

^ 48ttM|,' {> 
respectively. In this universal way the energy flows from the inflaton to ordinary particles. 
There is no any amplification of the energy drain to the SM particles due to coherent effects 
(like ones in Ref. [6]), since the produced particles interact strongly enough. 

Formulas (7), (8) generally mean that scalar particles play the major role in the scalaron 
decay process. The dominant contribution to the decay rate comes from coupling to kinetic 
term in Eq. (5). There is no similar term in the fermion case, see Eq. (6), because of the 
conformal invariance 3 of the fermion kinetic term. The energy transfers to the relativistic 
scalars by the time t reh > l/X^^V-w We define the reheating temperature T re ii as an 
effective temperature of produced relativistic matter at the time of equality between inflaton 
and matter energy densities. Then in the absence of coherent effects mentioned above, for N s 
scalar components contributing to the inflaton decay the effective temperature of produced 
relativistic matter is 

T reh * 4.5 x 10- 2 x . (^) 1/2 , (9) 

where g* denotes the number of relativistic species. In case of the Standard Model with 
= 106.75 and N s = 4, Eq. (9) gives numerically 

T reh « 3.1 x 10 9 GeV , (10) 

which is well in the region where the gauge interactions of the Standard Model are in equi- 
librium. So, the value (10) is the maximum temperature of the primordial plasma or the 
reheating temperature, indeed. 



2 These estimates are in agreement with similar ones in the Jordan frame originally obtained in [1, 5], 
while (7) and (8) are two times larger and four times smaller, respectively, than the estimates in Ref. [2]. 

3 This implies quite a similar situation for a scalar conformally coupled to gravity in the Jordan frame. 
We do not consider this case, but note that the coherent processes at preheating stage we alluded to above 
may be relevant then. 
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With low reheating temperature (10) the post-inflationary matter-dominated stage lasts 
long enough, so that the subhorizon scalar perturbations growing at this stage become non- 
linear (for a recent study see e.g. [7]). One expects that scalarons then start to form self- 
gravitating clumps of linear sizes much smaller than horizon, so the Universe still expands 
as at matter domination. The scalaron overdensity in clumps can be estimated similar 
to analysis of dark matter halos and clumps in the late-time Universe, see e.g. [8]. The 
scalaron density in clumps is not high enough to initiate scalaron scatterings due to self- 
interaction. Therefore, relativistic (SM) particles production can be described as decays of 
(non-relativistic) scalarons and the estimate of reheating temperature (10) remains intact. 

The decay rate formulas (7), (8) show that the scalaron produces ordinary particles 
in a universal way, so that all scalars and fermions of the theory will eventually populate 
the expanding Universe. With all non-gravitational interactions between ordinary particles 
switched off, the resulting abundances would be determined mostly by spin and, for fermions, 
by mass of the particles. 

Now let us consider a new field which is either free in the Jordan frame or couples to other 
fields very weakly, so that these new particles never equilibrate in the primordial plasma 4 . 
If stable at cosmological time scale, this particle is a good candidate to be dark matter. 

Given the discussion above, the particle mass is the only free parameter and hence its 
value is fixed by the requirement of comprising all dark matter whose relative contribution 
to the present energy density p c is [3] Q>dm ~ 0.223. The mass of dark matter particles mpM 
and their number density at present Udm,o are related to p c and Qdm as follows 



where we introduced the present entropy density s . 

Let us estimate the entropy-to-dark-matter ratio s/udm by the time of reheating. The 
dark matter production after inflation can be described as the decay of non-relativistic 
scalarons of mass whose number density evolves with scale factor a — a (t) at t > t re h 
as 



Here (f) re h, CL re h refer to the values of inflaton field and scale factor at reheating. We again 
neglect any coherent effects related to the new fields and also treat this new decay channel of 
scalaron as subdominant one, so the reheating temperature is still given by Eq. (10). Both 

4 In particular, the Planck-scale suppressed nonrenormalizable couplings are allowed. 



&DM Pc So 



(11) 



so nuM,o 
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points are discussed below in due course: we will see that the obtained results justify our 
choice. Assuming two-body decays of the inflaton to dark matter particles with decay rate 
T^dm one writes down the Boltzmann equation for the dark matter density, 

d 



— (n DM a 3 ) = 2n (a) T^ DM a 3 



This equation has a solution 



17>DM {Ueh) — ~ ^(p—tDM Ueh 



where p re h is the total energy density at the time of reheating, when the corresponding 
contributions of inflaton and relativistic species coincide. The energy density is related to 
the Hubble parameter at reheating H re h by 

p reh = 2g^T r \ h = 3M P H? eh . 

Taking for the numerical estimate of the reheating time 

1 

treh 



V3 H re h 

we finally obtain for the entropy-to-dark-matter ratio at reheating 



s (r re ,,) = ^#-^^. (12) 



n-DM 3 y/TE r^flM Mp 

Since this ratio remains intact at the further hot stages, we have the estimate 

s s 
noM noM,o 

Eqs. (12), (11) imply the following equation for the mass of dark matter particles universally 
produced by the scalaron decay, 

Qdm Pc 2tt yfgl p T reh 

m DM = - f— = -—- ■ 13) 

so 3V15 r^->DM M P 

Then in the cases of scalar (7) and fermion (8) one obtains the numerical estimates 
1.3xlO- 5 MpV / V^-V /2 / r 9* \ 1 IU^dm\[ P 



V P J W V106.75; VO.223^ V0-52xl0- 5 GeV/cm3y' ' 

(14) 

( P V /2 fN s \ 1/6 ( g* V/12 fn DM \ 1/3 f p c \ 

1.2x10 GeVx ^ 3xlQ _ 5M J i^j (— ) ^— j [ 52xl0 ^ GeV/cm3 ) 



(15) 
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Here we have corrected numerical coefficients (by tens of percent) to match the results ob- 
tained with numerical integration of the Boltzmann equations which consistently describe 
both the energy transfer from scalaron to all relativistic scalars and the dark matter produc- 
tion. 

The dark matter particles are produced highly relativistic with 3- momenta ~ fi/2, which 
exceed T re h and hence momenta of particles in the plasma by about four orders of magnitude, 

H ~ 10 4 x T reh . 

So, the free scalars of lOkeV mass would contribute to hot dark matter component only. 
This component is certainly subdominant, hence free scalars cannot solve the dark matter 
problem in the case of i? 2 -inflation. The viable choice of dark matter is free fermions of mass 
m DM ~ 10 7 GeV. These particles would form cold dark matter. 

We end up the discussion of universal dark matter in /^-inflationary model with several 
remarks. First, for the considered case of free boson particles coherent effects [6] during 
scalaron oscillations could change the obtained results. However, it only makes worse the 
situation with scalar dark matter, and we do not consider these processes. 

Second, if at a later stage of the Universe expansion the entropy gets produced (say, due 
to decays of nonrelativistic particles out of equilibrium, etc.), the estimates (14), (15) have 
to be corrected to account for the corresponding dilution factor r s = s new /s id- Thus, r.h.s. 
of (14) has to be divided by r s , while r.h.s. of (15) has to be divided by r^ 3 . In realistic 
models with moderate (if any) entropy production this does not save scalar dark matter, 
and only mildly changes the prediction for the fermion dark matter mass (15). 

Third, in the extensions of particle physics with additional scalars coupled to the SM 
fields, both reheating temperature (9) and would-be scalar dark matter mass increase. With 
very large number of new scalars this makes the scalar dark matter viable. 

Fourth, if there are other sources of dark matter in the model — other stable particles 
or other mechanisms of out-of-equilibrium production of heavy fermions — the estimate (15) 
implies the upper limit on the mass of universal dark matter in i? 2 -inflationary model. 

Fifth, in case of the Majorana fermions the mass is larger by a factor of 2 1 / 3 as compared 
to the Dirac case (15), since the decay rate of scalaron into the Majorana particles is two 
times smaller. 

Sixth, the estimates (14), (15) mean that heavier free particles are forbidden in models 
with i? 2 -inflation, otherwise the Universe would be overdosed. This conclusion is true for 
the particles lighter than the scalaron, when these estimates of particle production rates are 
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applicable. Clearly, production rate of heavier particles by means of scalaron oscillations is 
suppressed as compared to (7), (8). The heavier the particles, the stronger the suppression. 
Hence, free particles of mass m ^> fi are not forbidden in R 2 - inflation. Moreover, with 
strongly (most probably exponentially) tuned value of mass m ~ they can form viable 
dark matter. We do not study this situation here. 

3 Leptogenesis in i? 2 -inflation 

The study presented in the previous Section revealed that heavy free (in the Jordan frame) 
fermion is a viable dark matter candidate in i? 2 -inflationary model. These fermions are 
naturally produced in the post-inflationary Universe by scalaron decays, so the right amount 
of dark matter is achieved with the mass of about 10 7 GeV (15). One can further ask 
whether it is possible to make use of this universal production mechanism to unravel another 
cosmological problem which SM fails to solve: the baryon asymmetry of the Universe. The 
answer is positive, and we illustrate it with example of nonthermal leptogenesis via decays of 
heavy sterile neutrinos [9] universally produced by scalaron decays. This particular example 
is strongly motivated by observed oscillations of active neutrinos, phenomena which also lack 
an explanation within SM. Heavy sterile neutrinos give masses to active neutrinos via seesaw 
mechanism. This mechanism is responsible for the hierarchy between the neutrino masses 
and electroweak scale. In this Section we show that both the correct neutrino masses and 
successful leptogenesis can be realized in the i? 2 -inflationary model. 

The modification of SM we consider consists of two new Majorana fermions Nj, I = 1, 2, 
which are right singlets with respect to the SM gauge group. The most general renormalizable 
lagrangian for these fermions is 

Mr - 

C = iNa^Nj - y aI L a Nj$ - -^N^Nj + h.c. , (16) 

where y a i are new Yukawa couplings, $ is the SM Higgs doublet and $ = e $* with e being 
antisymmetric 2x2 matrix. 

When electroweak symmetry breaks down, $ acquires vacuum expectation value, $ T = 
(0,v/y/2), v = 246 GeV. Then (16) yields active-sterile mixing in the neutrino mass matrix. 
Assuming yv <C Mj, this mixing gives for the active neutrino masses: 

v 2 

m »af3 = -J2y«i^jfypi ■ ( 17 ) 
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Hence, m v <C v, that is the seesaw mechanism 5 . The formula (17) together with atmospheric 
neutrino mass splitting, Am 2 tot ~ 3 x 1CT 3 eV 2 [10], implies an order-of- magnitude upper 
limit on the lightest sterile neutrino mass: 

M Ughtest < 7 £== ~ 5 x 10 14 GeV . (18) 

The smaller the Yukawa couplings y a i and the heavier the active neutrinos, the smaller 
M UghtesU see Eq. (17). 

Lagrangian (16) added to SM explicitly violates the lepton number. The source of CP- 
violation, which is one of the Sakharov's conditions [11] for successful baryo- and lepto- 
genesis, is the complex Yukawa couplings y a j. The third Sakharov's condition, departure 
from thermal equilibrium for lepton number violating processes, is achieved due to decay of 
non-relativistic sterile neutrinos, which is evidently an out-of-equilibrium process. 

Let us discuss the production of these heavy sterile neutrinos in the post-inflationary 
Universe with oscillating scalaron. Neutrinos are produced via Yukawa-type coupling to 
scalaron, see Eq. (6). The production is effective only for neutrinos lighter than scalaron. 
Indeed, the amplitude of scalaron oscillation never exceeds M P . At the same time the 
scalaron coupling to a fermion is proportional to the fermion mass, so the latter remains 
almost intact. Therefore, an enhancement of heavy fermion production by Yukawa source 
due to periodic decrease of fermion effective mass [12] does not work here. 

Thus, we consider the models with sterile neutrinos lighter than scalaron, Mi < m^. By 
the time of reheating, t ~ t re h, the ratio of number density of produced in scalaron decays 
heavy neutrinos to entropy is given by Eq. (12). Substituting relevant numbers in (12) and 
(8) (recall that for Majorana fermions the decay rate is two times lower than for the Dirac 
fermions) one obtains 

^ (T reh ) = 2.9 x 1(T 6 x ( ^-—] 2 , (19) 

s K ' \5 x 10 12 GeVy ' V ; 

where we have included a correction of about tens of percent in order to match the solution 

of the relevant Boltzmann equations (see comment below Eqs. (14), (15)). 

These sterile neutrinos, however, decay quite rapidly due to the same Yukawa couplings 

in (16) responsible for the seesaw mechanism. Their total decay rates are 

r Ml I |2 

a 

5 Note that with only two seesaw sterile neutrinos we end up with one out of three active neutrino being 
massless. 
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Disregarding any hierarchy in the neutrino Yukawa sector one takes for the order-of-magnitude 
numerical estimate of Yukawa couplings relation (17) with m 2 ~ Am^ jm and finds for the 
lightest sterile seesaw neutrino 



jV 7 



An v 2 

Thus, sterile neutrino of mass above 10 10 GeV decays before reheating 6 . This is a very 
process which generates lepton asymmetry in the early Universe. To simplify formulas we 
assume further, that only the lightest neutrino, Aq, contributes and Mi <C M 2 . Then the 
value of lepton asymmetry A L = n L / s is 

A L = 5 L -^, (20) 

s 

with [9] (for reviews see e.g. [13]) 

SMi 1 ^ T / „ v 2 \ 3M 1 y/Am 2 atm 

= ~ w ' E^S7 \ V aimya2 ^ — ' ( } 

where we have used the order-of-magnitude estimate (17) with m 2 v ~ Am 2 tm . Finally, from 
(20), (21) and (19) we obtain 

A '~ h5 X 10 ' 9 X ( 5 x loS GeV ) 3 ' < 22 > 

The asymmetry is generated mostly at late times just before reheating, when the ratio 
of sterile neutrino (if stable) density to entropy becomes maximum, see Eq. (12). The 
sterile neutrinos are quite heavy so the lepton asymmetry is not washed out by inverse 
decay processes, which are out of equilibrium. Likewise, another important lepton violating 
process, scatterings $L -h- $Z through virtual sterile neutrinos (with cross section a ~ g^r) 
are out-of-equilibrium at the epoch of interest. Thus, lepton asymmetry (22) transfers later 
on into baryon asymmetry by sphaleron processes [14], 



A L /3 ~ 0.5 x 10~ 9 . 



Given the present baryon asymmetry of the Universe A^o = 0.88 x 10~ 10 [3] and the estimates 
above, we conclude that heavy seesaw sterile neutrinos of mass 10 12 -10 13 GeV, produced in 
scalaron decays, in i? 2 -inflationary model can be responsible for the baryon asymmetry of 



3 Eqs. (19) and (10) mean that sterile neutrino decays contribute not much to the Universe reheating. 
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the Universe. Since the sterile neutrino are lighter than scalaron, the available room is not 
large, that one could interpret as an order- of -magnitude prediction of the baryon asymmetry, 
which is otherwise (and usually) appears as an accidental number. 

The estimates above are an order-of-magnitude only. In the particular model under 
discussion with only two seesaw sterile neutrinos one can relate all the Yukawa couplings to 
the measured active neutrino mass differences and mixing angles, so that only few model 
parameters remain free. In this way one can refine our estimates. Likewise the numerical 
estimates may be changed in models with a particular hierarchy of Yukawa couplings in the 
neutrino sector. 

Note in passing that second, heavier sterile neutrino N2 if not too heavy can also con- 
tribute to lepton asymmetry generation, its decay products do not equilibrate. There is also 
some contribution from sterile neutrino decays at earlier times, when they are relativistic. 
The inverse processes are rare because the particle density is small. 

4 Discussion 

To conclude, we have shown that in i? 2 -inflation the scalaron decay responsible for reheating 
can also take care of dark matter production and generation of the baryon asymmetry of 
the Universe. The latter occurs through nonthermal leptogenesis by introducing seesaw 
neutrinos. Thus, neutrino oscillations are also explained in this model, which pretends to be 
phenomenologically self-contained. 

We end up with several comments. First, at relatively low energy scale, the model 
we discussed can be further modified in gravity sector (attempting to explain the late- 
time accelerated expansion of the Universe, as e.g. in Refs. [15, 16]) or in particle physics 
sector (say, by adding PQ-axion [17] to solve the strong CP-problem). All these and similar 
modifications can be safely adopted provided the early time cosmology remains intact. 

Second, we do not discuss here the gauge hierarchy problem in particle physics. Certainly, 
i? 2 -term does not contribute to this problem: all couplings of scalaron to other particles 
are suppressed by the Planck scale. With the absence of new scales in particle physics, 
one cannot discard that still unknown complete quantum theory at gravity scale (quantum 
gravity) is responsible for the cancellation of dangerous quantum corrections (if any) to the 
Higgs boson mass. In fact, a nontrivial space structure that shows up at Planck scales could 
also invalidate the strong CP-problem (see e.g. [18]). Thus, quantum gravity could solve 
all problems, indeed. Otherwise, one can think about other solutions usually considered. In 
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particular, our logic can be adopted for supersymmetric extensions of the SM and models 
with axion. In all these cases free fermion again serves as dark matter candidate, though it is 
somewhat heavier than in the case of SM because of larger number N s of scalars responsible 
for reheating and larger number of degrees of freedom g*, see Eq. (15). 

Third comment, which is related to the second one. The heavy sterile neutrinos (of mass 
Mat) we used contribute to gauge hierarchy problem making it worse as compared to the 
SM. Indeed, their coupling to the SM Higgs boson gives rise to both divergent and finite, of 
order Mjv, contributions to the Higgs boson mass. This worsening can be avoided if sterile 
neutrinos coupled to the Higgs boson have masses of order electroweak scale or smaller, as 
e.g. in z/MSM model [19]. Indeed, heavy sterile neutrinos can be replaced by the light ones 
in i? 2 -model we discuss. In that case, however, sterile neutrinos, required for leptogenesis, 
should be produced not by the scalaron decay (which, as we saw, requires large mass to be 
efficient enough), but via neutrino oscillations in the early Universe [20]. The i? 2 -model with 
one heavy fermion (free in the Jordan frame) and two light seesaw sterile neutrinos can be 
both phenomenologically and theoretically self-contained upto the quantum gravity effects. 
At the same time, this modification has also one more advantage as compared to what we 
discuss in the main text: sterile neutrino sector here can be directly tested [21] in particle 
physics experiments. 

We thank S. Odintsov, V. Rubakov, M. Shaposhnikov, A. Starobinsky, I. Tkachev and 
M. Vasiliev for valuable discussions. The work is supported in part by the grant of the Pres- 
ident of the Russian Federation NS-5525.2010.2 (government contract 02.740.11.0244). D.G. 
thanks ITPP EPFL and the organizers of the long-term workshop in Yukawa Institute YITP- 
T-10-01 for hospitality. The work of D.G. is supported in part by the Russian Foundation 
for Basic Research (grants 08-02-00473a), by FAE program (government contract n520) and 
by SCOPES program. The work of A. P. is supported in part by the "Dynasty" Foundation 
(awarded by the Scientific board of ICPFM), by FAE program (government contract n2598) 
and by the grant of the President of the Russian Federation MK-4317.2009.2. 

References 

[1] A. A. Starobinsky, Phys. Lett. B 91 (1980) 99; A. A. Starobinsky, "Nonsingular model 
of the Universe with the quantum-gravitational de Sitter stage and its observational 
consequences," in: Proc. of the Second Seminar "Quantum Theory of Gravity" (Moscow, 



12 



13-15 Oct. 1981), INR Press, Moscow, 1982, pp. 58-72 (reprinted in: Quantum Gravity, 
eds. M.A. Markov, P.C. West, Plenum Publ. Co., New York, 1984, pp. 103-128). 

[2] T. Faulkner, M. Tegmark, E. F. Bunn and Y. Mao, Phys. Rev. D 76 (2007) 063505 
[arXiv: astro-ph/06 12569] . 

[3] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180 (2009) 330 
[arXiv:0803.0547 [astro-ph]]. 

[4] S. Weinberg, The Quantum Theory of Fields , 3rd volume — Cambridge University 
Press, 1996. 

[5] A. Vilenkin, Phys. Rev. D 32 (1985) 2511. 

[6] L. Kofman, A. D. Linde and A. A. Starobinsky, Phys. Rev. Lett. 73 (1994) 3195 
[arXiv:hep-th/9405187]. 

[7] K. Jedamzik, M. Lemoine, J. Martin, JCAP 1009 (2010) 034. [arXiv: 1002.3039 [astro- 
ph.CO]]. 

[8] C. G. Lacey, S. Cole, Mon. Not. Roy. Astron. Soc. 262 (1993) 627-649; V. Berezinsky, 
V. Dokuchaev, Y. Eroshenko, Phys. Rev. D68 (2003) 103003 [astro-ph/0301551]. 

[9] M. Fukugita and T. Yanagida, Phys. Lett. B 174 (1986) 45; G. Lazarides and Q. Shah, 
Phys. Lett. B 258 (1991) 305. 

[10] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667 (2008) 1. 

[11] A. D. Sakharov, Pisma Zh. Eksp. Teor. Fiz. 5 (1967) 32 [JETP Lett. 5 (1967 
SOPUA,34,392-393.1991 UFNAA,161,61-64.1991) 24]. 

[12] G. F. Giudice, M. Peloso, A. Riotto and I. Tkachev, JHEP 9908 (1999) 014 [arXiv:hep- 
ph/9905242]. 

[13] W. Buchmuller, R. D. Peccei and T. Yanagida, Ann. Rev. Nucl. Part. Sci. 55 (2005) 
311 [arXiv:hep-ph/0502169]; S. Davidson, E. Nardi and Y. Nir, Phys. Rept. 466 (2008) 
105 [arXiv:0802.2962 [hep-ph]]. 

[14] S. Y. Khlebnikov and M. E. Shaposhnikov, Phys. Lett. B 387 (1996) 817 [arXiv:hep- 
ph/9607386]. 



13 



[15] S. Nojiri and S. D. Odintsov, Phys. Rev. D 68, 123512 (2003) [arXiv:hep-th/0307288]; 
S. Nojiri and S. D. Odintsov, eConf C0602061, 06 (2006) [Int. J. Geom. Meth. Mod. 
Phys. 4, 115 (2007)] [arXiv:hep-th/0601213]; S. Nojiri and S. D. Odintsov, "Non-singular 
modified gravity unifying inflation with late-time acceleration and universality of viscous 
ratio bound in F(R) theory," arXiv: 1008.4275 [hep-th]. 

[16] A. A. Starobinsky, JETP Lett. 86 (2007) 157 [arXiv:0706.2041 [astro-ph]]. 

[17] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38 (1977) 1440. 

[18] S. Y. Khlebnikov and M. E. Shaposhnikov, Phys. Lett. B 203 (1988) 121; S. Khlebnikov 
and M. Shaposhnikov, Phys. Rev. D 71 (2005) 104024 [arXiv:hep-th/0412306]. 

[19] T. Asaka, S. Blanchet and M. Shaposhnikov, Phys. Lett. B 631 (2005) 151 [arXiv:hep- 
ph/0503065]; T. Asaka and M. Shaposhnikov, Phys. Lett. B 620 (2005) 17 [arXiv:hep- 
ph/0505013]. 

[20] E. K. Akhmedov, V. A. Rubakov and A. Y. Smirnov, Phys. Rev. Lett. 81 (1998) 1359 
[arXiv:hep-ph/9803255]. 

[21] D. Gorbunov and M. Shaposhnikov, JHEP 0710 (2007) 015 [arXiv: 0705. 1729 [hep-ph]]. 



14 



